Quantum statistical properties of the radiation field in a cavity with a movable mirror 
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A quantum system composed of a cavity radiation field interacting with a movable mirror is 
considered and quantum statistical properties of the field are studied. Such a system can serve in 
principle as an idealized meter for detection of a weak classical force coupled to the mirror which 
is modelled by a quantum harmonic oscillator. It is shown that the standard quantum limit on 
the measurement of the mirror position arises naturally from the properties of the system during 
its dynamical evolution. However, the force detection sensitivity of the system falls short of the 
corresponding standard quantum limit. We also study the effect of the nonlinear interaction between 
the moving mirror and the radiation pressure on the quadrature fluctuations of the initially coherent 
cavity field. 
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I. INTRODUCTION 

Due to recent experimental efforts to detect gravita- 
tional waves the problem of fundamental limitations 
on the detection sensitivity imposed by the quantum na- 
ture of the measurement device has become of practical 
importance. The research during the last decade focused 
mostly on interferometric schemes in which the relative 
shift in the positions of two end mirrors, caused by a 
gravitational wave, results in a phase shift between the 
light beams in the two arms. This phase shift can be de- 
tected by measuring the output light. The first analysis 
of quantum limitations on the measurement sensitivity in 
such an interferometric scheme was made by Caves |2j in 
the early eighties. According to this analysis, there are 
two sources of quantum noise: the photon-counting noise 
and the radiation-pressure noise. When these two sources 
of noise are balanced by adjusting the light intensity, the 
sensitivity achieves the so-called standard quantum limit 
(SQL). 

More recently, Pace et al. |3) studied quantum limita- 
tions on the interferometer sensitivity using a more de- 
tailed analysis based on a Hamiltonian description of the 
interaction of light with the movable end mirrors. Some 
other works [f|-{7| were devoted to the study of quantum 
noise reduction in simpler optical systems: a light beam 
reflected by a movable mirror or a single Fabry-Perot cav- 
ity with one movable mirror. Assuming high light intensi- 
ties, the above works mostly employed semi-classical cal- 
culational techniques, with the quantum fluctuations of 
the light field considered as a small addition to the large 
classical amplitude. The quantum theory of the radiation- 
pressure fluctuations on a mirror was also studied in Q . 

In the present paper we consider a simple but neverthe- 
less instructive system, consisting of a cavity radiation 
field interacting with a movable mirror |9|dTlf|. As was 
shown recently |l^[L3|, such a system can be utilized to 
generate a variety of non-classical states (especially, the 
so-called Schrodinger-cat states) of both the cavity field 
and the mirror. We employ the fully quantum description 



of the system and study quantum statistical properties of 
the cavity field. While the photon statistics of the field 
does not change during the interaction, the phase dis- 
tribution and quadrature fluctuations exhibit interesting 
dynamics. 

If the mirror (modelled by a quantum harmonic oscilla- 
tor) is driven by an external classical force, this will result 
in a phase shift of the cavity mode. We study quantum 
limitations imposed on the detection of the force by the 
intrinsic quantum uncertainty in a measurement of this 
phase shift. It is satisfactory to find that the SQL for 
the measurement of the mirror position naturally emerges 
from the dynamical properties of the system. On the 
other hand, the sensitivity of the force detection in this 
scheme does not reach the corresponding SQL. We also 
show that analytical approximations made to describe 
phase properties of the field should work very well for 
light intensities needed to achieve the SQL in appropriate 
experimental setups. 



II. THE STANDARD QUANTUM LIMIT IN THE 
SIMPLEST OPTOMECHANICAL SENSOR 



As is well known, the SQL arises due to the Heisenberg 
uncertainty relation for the position and momentum of 
the movable mirror. This limit can be understood con- 
ceptually by considering a simple measurement scheme 
consisting of a movable mirror (modelled, for example, by 
a linear harmonic oscillator of mass m and angular fre- 
quency u) m ) and a monochromatic light beam of angular 
frequency c^o- A shift z in the mirror position, caused by 
an external classical force (e.g., by a gravitational wave), 
results in a phase shift <fi = 2ujqz/c of the reflected light 
beam (here c is the velocity of light). Such a scheme is 
referred to as an optomechanical sensor. The experimen- 
tal progress towards the observation of quantum limits 
in such a sensor was reported recently in j^J. For a re- 
cent theoretical treatment of such a scheme, including a 
detailed analysis of various sources of noise, see |15) . 
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If the incident light beam is in a coherent state with 
mean photon number N, the phase uncertainty is A</> ~ 
l/(2iV 1 / 2 ) (provided that N > 1). This phase-fluctuation 
noise results in the uncertainty (Az) p f ~ c/(Au)qN 1 ' 2 ) 
in the measurement of the mirror position. Increasing 
the light intensity, one can suppress the phase-fluctuation 
noise, but this leads to increased fluctuations in the 
radiation pressure which on the average transfers the 
momentum (2hujo/c)N to the mirror. These radiation- 
pressure fluctuations result in the mirror position uncer- 
tainty (Az) rp ~ 2HujqN 1 ^ 2 /(mcu) m ). Optimizing the sum 
(Az) 2 = (Az)p f + (Az) 2 p as a function of N, one finds 

the uncertainty (A^) opt which is y/2 times the SQL for a 
harmonic oscillator, 



(Az) S q L = ^h/(2maj m ). 
The optimum mean number of photons is 
N opt = mw m c 2 /(8fiu 2 ). 



(1) 



(2) 



In many applications (e.g., for detection of gravitational 
waves) it makes sense to consider what limitations apply 
due to the quantum nature of the measuring apparatus 
on the detection of a weak classical force coupled to the 
oscillator pq |. Consider a constant force F that acts on 
the mirror during time t which is much shorter than the 
oscillator period 2ir/uj m . In order to detect this force, the 
momentum pp — Ft it transfers to the mirror should be 
larger than the SQL for the measurement of the oscillator 
momentum, (Ap)sql = (fimWm/2) 1 / 2 . Correspondingly 
the minimum detectable force is 



-Fsql = t 1 (Ap)sql = t 1 \Jhmu} m j2. 



(3) 



It was argued Jl7|-p0| that the SQL on the detection of 
a weak classical force (which causes the displacement of 
the movable mirror in the sensor) can be surpassed using 
correlations between the phase-fluctuation noise and the 
radiation-pressure noise. In particular, it was predicted 
[Q-Q that amplitude-phase correlations, created by the 
nonlinear interaction between the movable mirror and the 
radiation pressure, will produce squeezing in the reflected 
beam. Consequently, it was proposed Q to use homodyne 
detection of the squeezed quadrature in order to compen- 
sate for the radiation-pressure noise, thereby making it 
possible to overcome the SQL with high enough light in- 
tensity. 



III. THE MODEL 

The solution of the wave equation for the radiation 
field in a cavity with moving boundaries is in general 
a very complicated problem even in the classical case 
P,^0|. However, in the adiabatic approximation, when 
the amplitude and frequency of the mirror motion are 
much smaller than the wavelength and frequency of the 



cavity mode, the dynamics of the system can be described 
by an effective Hamiltonian jy]]. Of course, all phenom- 
ena associated with the dynamical Casimir effect (e.g., 
the resonant generation of photons) are neglected in this 
description. 

The free Hamiltonian for the system of the radiation 
field in a linear cavity with one movable mirror is 



Hq = H c + H rn = hio c a) a + fkj m b*b, 



(4) 



where u c and a are the angular frequency and the boson 
annihilation operator of the cavity mode, while u> m and 
b are the angular frequency and the annihilation opera- 
tor of the movable mirror modelled by a quantum har- 
monic oscillator. For a cavity of free length L, one has 
lo c = n(c/L)n, n £ N. If the mirror displacement from the 
equilibrium position is z, then the cavity length changes 
L — > L + z, and the cavity angular frequency changes as 



'"I 



(5) 



Therefore, the effective interaction Hamiltonian reads 



-huj c —a^ a = —hga) a(b + b^), (6) 



where 



/ h \ 1/2 
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\2mujr, 

/ * \ V 2 
( 

L \ 2mu> r , 



(7) 
(8) 



Here, m is the mass of the movable mirror. If, in addition, 
we assume that the mirror is driven by a classical external 
time-dependent force F(t), the corresponding part of the 
Hamiltonian is given by 



H F = -F(t)z = -hf(t)(b + tf), 
f(t) = {2muo m h)- 1 ' 2 F{t). 



(9) 
(10) 



Then the total Hamiltonian H to t = Hq + H[ nt + Hp reads 
H tot = nw c a)a + fiw m bh - %a f a(6 + 6 f ) - hf(t)(b + tf). 

(11) 

It is known that the time evolution operator for this 
system can be written in a closed form |l^ , [l3|] . One can 
easily verify that this evolution operator is given by 

U (t) = cxp(— ira/ar) exp[— ib^br + i(ka ! a + \)(b + b ')t]. 

(12) 

With some algebra jl3| , this operator can be rewritten as 

U(t) = exp(— ira^ar) exp[j/i(fca^a + A) 2 ] 

x cxp[(fca t a + A)(7 7 6 t - rfb)] exp(-i6 t 6r). (13) 
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Here, r = oj m t is the scaled time and other scaled param- 
eters are r — u> c /u> m and k = g/u> m . Also, A, /i, and r\ are 
dimensionless functions of time: 



A(r) 



1 /""V" 

T 

fi(r) = t — sin r, 
?7(t) = 1 — exp(— ir) 



f(t')dt', 



(14) 

(15) 
(16) 



An interesting feature of this system is that for times 
t = 2ir, 4-7T, . . . (i.e., after a full cycle of the mirror oscil- 
latory motion), one finds rj = 0, and the field and mirror 
subsystems become fully disentangled (if they were mi- 
correlated initially). Furthermore, at these moments, the 
mirror returns to its initial state. This effect is quite gen- 
eral as it does not depend on the initial states of the sub- 
systems. Another interesting property is that the photon 
statistics of the cavity field does not change during the 
interaction. 

The existence of the closed form ( |l2| ) for the evolution 
operator can be understood if we note that the Hamil- 
tonian (HJ) is given by a linear combination of operators 
which close an algebra. First, recall that a solvable Lie 
algebra (the so-called oscillator algebra) is closed by the 
four operators {N b , b, b\ I}, where N b = b^b is the num- 
ber operator and I is the identity operator. Non- vanishing 
commutation relations are 

[N b ,b} = ~b, [N b ,tf] = tf, [b,tf]=I. (17) 

Alternatively, one can choose for the oscillator algebra 
the fully Hermitian basis {N b , X b , Y b , I}, where X b = 
2~ 1 / 2 (b + &t) and Y b = 2-^ 2 {b - are the scaled 

position and momentum operators. Non- vanishing com- 
mutation relations are 



[N b ,X b ] = -iY b , [N b ,Y b ]=iX b 



[X b) Y b ] = il. 

(18) 



It is straightforward to see that the oscillator algebra can 
also be constructed using two bosonic modes a and b. The 
Hermitian basis is given in the two-mode realization by 
{N b , N a X b , N a Y bl N 2 }, where N a — a^a is the number 
operator for the mode a, and non- vanishing commutation 
relations arc 



[N b ,N a X b ] = -iN a Y b , 
[N a X b ,N a Y b ] =iN 2 a . 



[N b ,N a Y b ]=tN a X b , 



(19) 



More generally, it is also possible to consider a solv- 
able Lie algebra which is closed by the eight operators 
{N b ,N a X b ,N a Y b ,X b ,Y b ,N%,N a ,I} with non-vanishing 
commutation relations given by the combination of (|l8|), 
(H) and 



[N a X b ,Y b ] = [X b ,N a Y b ]=iN a . 



(20) 



Since the Hamiltonian ( pT[ ) is given by a linear combina- 
tion of the operators N a , N b , N a X b , and X b , the Baker- 
Campbell-Hausdorff theorem |2l|] assures the existence of 
a closed expression for the time evolution operator. 



IV. THE EVOLUTION OF COHERENT STATES 

Let us consider a situation in which both the cavity field 
and the mirror oscillator are initially prepared in coherent 
states with coherent amplitudes a and /?, respectively: 



l*(0)> = \a) 



a' 1 



n=0 



n) c (g) 



(21) 



In fact, a quite reasonable choice is to assume that the 
mirror is initially in the vacuum state (i.e., (3 = 0), but 
first we will consider an arbitrary [3 for the sake of gener- 
ality. Acting on the initial state ( pl| ) with the evolution 
operator (13), one obtains: 



|*(r)) 



— 2TT \ n 



71—0 V ^ 

3 (fcn+A)(7)6 t -?7*6) 



i^- 



Utilizing properties of the displacement operators p2[ , 
the state of the system may be written as (we omit an 
unimportant overall phase factor) 



l«| 2 /2 



y-\ [ae 

n=0 



_ *C( T )l n 



xe^ k2n2 \n) c ® | 7 „(r)) m . 



(23) 



In the reference frame rotating with the field frequency 
u/ c , the phase ^(r) is given by 

C(r) = 2fcA(r)^(r) + fclm \J3t](t)]. (24) 

In Eq. (|23| ), |7«(t)) is a coherent state with the amplitude 

ln (r)=0e- iT +[kn + X(T)]r,(T). (25) 

For t = 2-7T, 4-7T, . . ., the dependence of this amplitude on 
n disappears, and the system returns to an uncorrelated 
product state. 

In what follows we will be interested in quantum statis- 
tical properties of the cavity field. These properties can 
be calculated from the reduced density matrix, obtained 
by tracing out the mirror degrees of freedom, 

p c (r)=Tr m {|*(r))(*(r)|} 

= e -H 2 y-V- (ae<)"(a*e-<)"' 

n=0n'=0 



V n\n'\ 

xe^ fe2 ^- nra )(7„'|7n)|n)cc<n'|, 

where 

(ln'\7n) = exp (7*,7„ - i|7„| 2 - \\"fn>\ 2 ) 

= cxp[-/ifc 2 (n-n') 2 ] for (3 = 0. 



(26) 

(27) 
(28) 



Note that fi = ^\rj\ 2 = 1 — cost. Obviously, the Pois- 
sonian photon statistics of the cavity field is preserved 
during the evolution. 
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V. PHASE PROPERTIES OF THE CAVITY FIELD 

As can be seen from equations ( p3| ) and (|24|) , the exter- 
nal force acting on the mirror produces a phase shift of the 
cavity mode. At least in principle, this phase shift can be 
measured, thereby revealing information about the classi- 
cal force. Accuracy of such a measurement will inevitably 
depend on the phase properties of the radiation mode, be- 
cause an externally induced phase shift is detectable only 
if it is larger than the intrinsic phase uncertainty of the 
field. 



A. Canonical and realistic phase distributions 

In quantum mechanics, one can describe a measure- 
ment by the corresponding positive operator-valued mea- 
sure (POVM) ||§. The so-called canonical POVM for an 
idealized phase measurement is given by f24|j25|| 



IL 



|n). 



(29) 



n=0 



The corresponding phase distribution for the cavity field 
is 

Pca„(0) = nrpW^pJ = {2^Y 1 {9\ Pc \e). (30) 
Substituting expression ( |26| ) for p c into (^0|), we obtain 



|2 OO OC 



2tt 



EE 



n+n' 



— £ — <p Q )(ro— n') 



n=0 n'=0 
Jnk 2 (n 2 — n' 2 ) 



\Jn\n'\ 

{ln'\ln), 



(31) 



where ip a = arg(a). For j3 — 0, we can use expression 
( p8| ) for (7„/ 17„). Then the double summation in Eq. ( fjl| ) 
can be rearranged in such a way that P ca ,n{d) takes the 
form of a Fourier series: 



Pcan(0) 



1 OC 

- y 



e- 1 * q A q e l 



where 



An 



q— — CG 



N-4,1 £ 



c2n 



' Q y/n\{n+\q\)\' 



£ q = \a\ cxp (-ipk 2 q). 



(32) 

(33) 
(34) 



The problem of the quantum description of opti- 
cal phase and the problem of its measurement were 
widely discussed during the last decade (see, e.g., reviews 
p6| , p7fl ). However, existence of an experimental procedure 
for phase measurement with the canonical POVM ( |29| ) is 
still an open question. On the other hand, a possible way 
to measure the phase of the radiation field mode is by 
means of heterodyne detection (see, e.g., p8[). Relating 



to our model, a possible experimental procedure might 
be as follows. At some moment t one stops the interac- 
tion between the field and the movable mirror (by fixing 
the mirror position) and lets the radiation leak from the 
cavity. The output field is detected using the heterodyne 
scheme, which employs a strong local oscillator highly de- 
tuned from the signal. The two Fourier components of the 
photocurrent are proportional to the quadratures of the 
cavity field at the moment t, which can be used to deter- 
mine the phase. The POVM for this phase measurement 



n Q (0) = i [°°\A)(A\rdr, A 
k Jo 



JO 



(35) 



where \A) is a coherent state with amplitude A. The cor- 
responding phase distribution is obtained by integrating 
the Q distribution function, 



Q(A) = ~(A\p c \A) 

IT 

over the radial coordinate in the complex plane: 



(36) 



p OC 

P Q {6) = T r [U Q (9)p c ] = / Q{re ld )rdr. (37) 
Jo 



For the cavity field mode with p c of equation fl26|), we 
obtain 



2 OO OO 



n + n' + 2\ \a\ n+n ' 



xe - J (e-c-^)(«-n') e vfe 2 (« 2 -n' 2 )^ n/ | 7n ^ ( 38 ) 

We see that the distributions P ca n(#) and Pq{Q) differ by 
the factor 



1 f n + n' + 2 



\Jn\n'\ 



(39) 



in the double sum over n and n' . Correspondingly, Pq(0) 
is usually slightly wider than P C!al {9). (It has been argued 
|29| l that this widening is related to contamination by ex- 
ternal noise during the heterodyne detection.) For (3 = 0, 
the rearrangement of the double summation in ( |38|) gives 

-. 00 

P Q (0) = 5- E e-^^B^ 8 -^^, (40) 



where 



« ...-i.iyi«iy r (n + j|g| + i) ff' 

q ' Kq ho nn+\q\ + l) n\ ■ 



(41) 



and £ q is given by Eq. (p4|). Note that Bo = 1- The sum- 
mation in (|4l]) is proportional to the Kummer series [ |30[ , 
so we can write 



B a = 



i r (3M + i)^i 



q r(M + i) 



<&(iM + i,M + i ; a. (42) 
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Here, $(a, b; x) is the Humbert symbol for the confluent 
hypergeometric function (the Kummer function). Using 
properties of special functions, B q can also be written as 



'IT 



B q = ^-exp(-|a| ')«p(£/2)£, 



2 

x [h\i 



l)/2(^g/2) +I{\ q \+l)/2 



(e 9 2 /2)] 



(43) 



where /^(x) is the modified Bessel function of the first 
kind. In the strong-field limit \a\ 2 3> 1, we can use the 
asymptotic expansion for the confluent hypergeometric 
function or for the modified Bessel function to obtain (for 
fixed q): 



B a 



exp( 



(44) 



For high field intensities, the behavior of both phase dis- 
tributions is qualitatively very similar. In what follows 
we will consider the Pq(0) distribution, because of two 
reasons: (i) phase measurements with the corresponding 
POVM Hq(9) are experimentally feasible and (ii) the ex- 
istence of a closed expression for B q gives a calculational 
advantage. 

Let us assume the phase range —ir<0<ir. Then 
moments of the phase distribution are determined by 



0P Q (0)d0, 2 
A simple calculation gives 
! Im 



2 P Q {0)d0. (45) 



2 £ f 



oo 



1 



q=l H 



Re 



,-»(C+¥>a-Hr)g 



(46) 



(47) 



for (x(t) and (i(t). The leading term in /i(r) is r 2 /2, so 



the terms in the series (40) are weighted by the prefac- 
tor exp(— k 2 T 2 q 2 /2). Correspondingly, the contribution of 
terms with q much larger than q m = (fcr) -1 is rather in- 
significant. Then, for high field intensities (i.e., |a| much 
larger than q m and \a\ 2 S> 1), we may approximate B q 
by just the leading term in its asymptotic expansion (44). 
With all these assumptions, we derive 



Pq{0) 



1 

27 



E 

q— — oo 



exp[-io-V]exp[i(0-%], (48) 



a = kry/1 + £i, 



= C + Pa + £2; 



9 



T 

12 



T 

360' 



£2 = 



(49) 



(50) 



We see that the phase distribution Pq (8) is approximated 
by a Fourier series with Gaussian coefficients, 



exp[ 



1 2 2l 



P Q (8)e- l{e - §)q d0. 



(51) 



Though we did not find a closed expression for the series 
(|4q), a further approximation may be made. It is clear 
that the width of the distribution is of the order of a and 
it is centered near 9. When the distribution is very narrow 
(i.e., a <§; 1) and its center is far from the boundaries of 
the interval [—it, tt) (in comparison with the distribution 
width a), it is possible to extend the integration range 
m © from — oo to +oo without changing the result sig- 
nificantly. This is equivalent to the replacement of the 
Fourier series in (|4^) by the Fourier integral. Formally, 
we can apply the Poisson summation formula (see, e.g., 
§1) to the series <M) to obtain 



P Q (0) » (2w 2 ) 



-1/2 



exp 



2mirf/a 2 



Using these expressions, we can calculate the phase un- 
certainty, A0 = (0 2 — 2 ) 1 / 2 . This can be done numer- 
ically, but some important information may also be de- 
duced from an analytical approximation. 



B. Analytical approximations and numerical results 



(52) 

For a < 1, the series ([32]) converges much faster than 
d4g|). If only the leading term in (52) is left, the phase 
distribution is approximated by just the Gaussian, 



P Q {0) w (2w 2 )- 1/2 exp [-1(0 - 0) 2 /a 



(53) 



In our considerations we neglected losses which occur 
in a realistic case due to the relaxation of the cavity field 
and the dissipation of the mirror oscillator. For a cavity 
with a very high finesse and a mirror oscillator with a very 
high quality factor, neglecting the losses can be a fair ap- 
proximation for short interaction times, say, for t < 10~ 3 
s. Therefore, for low-frequency oscillators (with u> m ~ 2-7T 
Hz, as in the LIGO intcrferometric gravitational-wave de- 
tectors Q), it is reasonable to consider r < 10~ 2 . For 
such small values of r we can use approximate expressions 



centered about 6 and with the dispersion a. 

The accuracy of our approximations may be verified by 
calculating numerically the mean phase and the phase 
uncertainty A0 from the distributions ( |4C| ) and (48). For 
simplicity, we took £ = ip a — in these calculations. 
Figure 1 shows how the phase uncertainty A0 depends 
on the intensity of the cavity field. The approximation 
pq ) works extremely well for sufficiently large field in- 
tensities. The Gaussian approximatio n (J53|) is also excel- 
lent. In principle, A0 obtained from (|48[) coincides with 
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a except for the neighborhood of points where the mean 
phase is ±7r mod 2tt. The corresponding deterioration 
of the phase sensitivity is an artifact of the definition of 
the phase range from —tt to tt. The dependence of the 
mean phase 9 on \a\ is shown in Fig. 2. It is seen that 
the value 9 obtained in the Gaussian approximation dif- 
fers very insignificantly from the exact value of 9. One 
can also verify that the peaks in Fig. 1 appear exactly at 
the points for which 9 in Fig. 2 approaches tt. Figure 3 
shows the dependence of the phase uncertainty A9 on the 
scaled measurement time r. Once again, the approxima- 
tions work very well. 

We should note that the value of k we used (k = 7.0) 
has no special physical meaning. In fact, this value is 
much larger than one would expect in an experiment. We 
chose such a large k only in order to deal with smaller 
values of |a| and thereby avoid overflows in numerical 
calculations. However, this artificial scaling does not dis- 
tort the physical behavior of the system. Let us con- 
sider what will happen for a realistic value of k (e.g., 
t~6x 10~ 8 for the LIGO detectors). Then one obtains a 
much larger value of q m = (kr)^ 1 (e.g., g ra ~3x 10 9 for 
the LIGO detectors with the cavity storage time about 
10~ 3 s). Therefore, our approximations might be used, 
but for very large values of \a\. In fact, a simple calcu- 
lation gives = (4L/ci) 2 7V opt , where t is the measure- 
ment time and ^V pt is the mean photon number needed 
to achieve the SQL in an optomechanical sensor, as given 
by Eq. (||). With the LIGO parameters, one obtains 
-^Vopt ~ 2.5 x 10 21 , so our approximations will work per- 
fectly for the mean photon numbers N — \a\ 2 of the order 
of iVopt an d even smaller. 



VI. QUANTUM LIMITATIONS 

The principal result is that, in the region of small r, the 
phase uncertainty A9 is limited from below by the value 
kr (in fact, this limit is slightly higher due to the factor 
\/l + £\ in er, but this small correction is not important 
for understanding the basic behavior of the system). Now, 
let us consider what this phase uncertainty means regard- 
ing the measurement of the mirror shift. A phase shift 4> 
produced due to the mirror position shift z is given by 
4> ~ (2Bu> c /c)z, where 



B 



tc 



(54) 



is the effective number of bounces by light on the movable 
mirror during the measurement time t. If the uncertainty 
of the phase shift measurement is A9, the corresponding 
uncertainty of the mirror position shift is 

- L ■ (55) 



Az 



2Buj r 



-A9 



-A9. 



Since in our model the phase uncertainty A9 is limited 
from below by kr, the corresponding limit on the posi- 
tion shift uncertainty is 



(A*)r 



w m i gL 
kr = — = 



2mu) Tl 



(56) 



which is exactly the SQL for the position measurement of 
a harmonic oscillator. 

The limitation on the detection of an external classical 
force F(t) will depend, of course, not only on the uncer- 
tainty A9 of the phase measurement, but also on the value 
of the phase shift = 2/cA/i produced by this force. The 
corresponding signal-to-noise ratio is 



S _ ( f (t) _ 2kX(T)fi(r) 



N A9 A0 
For t -C 1, taking for A9 its lower limit kr, we obtain 



(57) 



S fcA(r)r 3 /3 



N 



kr 



(2hmuj m y 1/2 



"A 



3 J 



F(t')dt'. 



(58) 



Clearly, the sensitivity of the force detection depends on 
the total phase shift produced by the force during the 
time t = T/u> m . This result just reflects the fact that in 
our model the information about the system is inferred 
from a "single-shot" phase measurement at the time t, 
and not from a continuous monitoring of the system. In 
the case of a force F which is constant during the time t, 
we obtain 



«? , Ft 2 

A ~ (2hmu m )-^^- 
N 3u m 



(59) 



The minimum detectable force -Fmm is determined from 
the condition S/N = 1. Then we obtain 



18hm 



t 4 



6 

-Fsql- 



(60) 



So, with the method presented here, the minimum de- 
tectable force is quite far from the SQL. This result shows 
that the genuine sensitivity in a specific scheme can dif- 
fer significantly from the commonly accepted SQL for an 
abstract oscillator. 



VII. QUANTUM FLUCTUATIONS OF THE 
FIELD QUADRATURES 



It is well known (see, e.g., 32 ) that homodyne detec- 
tion with a strong local oscillator can be used to measure 
the rotated field quadrature 



(61) 



The angle ip is the phase of the local oscillator. A quan- 
tum state of the signal field a exhibits the phenomenon 
of squeezing if the uncertainty AX V = {(X 2 ) - (X^) 2 ) 1 / 2 
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is below its vacuum value (AY v ) vac 
ip. Since 



1 for some angle 



\ + 2a)a + a 2 e- 2iip + a) 2 e 2lip , 



(62) 



the uncertainty AX V at any time t can be determined by 
calculating the expectation values of a and a 2 (recall that 
the mean number of photons (a^a) = \a\ 2 is constant). 
For /3 = 0, a straightforward calculation gives: 

\a\ exp[i(C + (fi a + ^k 2 )} 



(a) = (a*)' - 



x exp 



-fik 2 - |a| 2 (l 



{a 2 ) = (a t2 )* = \a\ 2 exp[2i(C + <p a + 2/zfc 2 



x exp 



-4/i/c 2 - |a| 2 (l 



(63) 



(64) 



Using these results, we can study numerically the behav- 
ior of the quadrature uncertainty. Also, in the relevant 
limits \a\ 2 3> 1 and r -C 1, an approximate expression is 
obtained 

(AX V ) 2 « 1 + 2\a\ 2 (1 + e" 2r cos 2tf - 2e" r cos 2 i?) , 

(65) 

(66) 
(67) 



r = fc 2 T 2 (i + fc 2 H 2 r 4 /9), 
i? = C + ^ + fc 2 M 2 T 3 /3- 



The variance in Eq. ( p5[ ) is minimized for = n7r (n 
0,±1,±2,...). Then we find 



(AX V ) 2 «1 + 2H 2 (1- 



,-r\ 



(68) 



which is always larger than 1. 

Figures 4 and 5 show the basic features of the behav- 
ior of the quadrature uncertainty. For simplicity, we took 
C — Va = m the numerical calculations. In Fig. 4 the 
dependence of AX V on the homodyne phase tp is shown 
for various field intensities (the range of ip from zero to 
7r is chosen because AX V is a periodic function of ip with 
the period of n). For the parameters taken, the approxi- 
mation ( |65| ) is excellent, and the quadrature uncertainty 
AX V reaches a minimum for ip about fc 2 |a| 2 r 3 /3 (i.e., for 
$ = 0). The value of AX V at this minimum is plotted in 
Fig. 5 versus |a| for several interaction times. Some dis- 
crepancy between exact and approximate results appears 
for small values of r, but the approximation ([38|) is excel- 
lent for field intensities large enough to assure T > 1. For 
the parameters we considered, no squeezing of the field 
quadrature was found. 



paper might be immediately used as a practical sensor of 
weak classical forces (e.g., for detection of gravitational 
waves). Our study is aimed at understanding the most 
basic properties of the quantum dynamics of the radiation 
field in a cavity with a movable mirror. The evolution of 
the field phase distribution, studied in the framework of 
a fully quantum approach, is of interest because it helps 
to clarify how the information about an external force is 
carried by the light field in much more complicated realis- 
tic measurement schemes. In fact, it is very instructive to 
see how the evolution of the field phase uncertainty deter- 
mines the sensitivity of the scheme and naturally implies 
the SQL on the measurement of the mirror position. 

It should be emphasized that in our model we assumed 
that the information about the phase shift accumulated 
by the cavity field is read by means of a "single-shot" 
phase measurement. Of course, this only gives informa- 
tion about the integrated effect of an external force during 
the measurement time. However, if one wants to deter- 
mine the time dependence of the force, this will require a 
continuous monitoring of the field phase. In a consistent 
quantum description, one should take into account that 
the evolution of the system under the continuous mea- 
surement will be seriously affected by the measurement- 
induced state reduction. A model of continuous broad- 
band measurement for monitoring the position of a free 
mass was recently analyzed in p3[ . It would be also in- 
teresting to have such a kind of quantum analysis for the 
continuous measurement on the system of the radiation 
field interacting with a movable mirror. 
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VIII. DISCUSSION 

It should be clearly understood that we have no inten- 
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FIG. 1. The phase uncertainty A9 versus \a\ for k — 7.0 
and t = 0.01. The exact result is shown with the solid line 
and the approximate result, corresponding to the distribution 
is the dashed line. The difference between exact and 
approximate results appears only for \a\ < 300. If one ignores 
peaks, the corresponding curve coincides perfectly with the 
Gaussian dispersion a. 
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FIG. 2. The mean phase 8 versus |q| for k = 7.0 and 
r = 0.01. The exact result is shown with a solid line and 
the Gaussian approximation 9 — £2 with a dashed line. The 
difference between exact and approximate results is very in- 
significant. The mean phase 9 obtained from the approximate 
distribution (0) is not shown here because it is almost indis- 
tinguishable from the exact value. 
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FIG. 3. The phase uncertainty A8 versus the scaled time 
r for k = 7.0 and \a\ = 500. The exact uncertainty and 
the approximate one, corresponding to the distribution (pjj|), 
are indistinguishable. If one ignores peaks, the corresponding 
curve coincides perfectly with the Gaussian dispersion a. 




FIG. 4. The quadrature uncertainty AX V versus <fi for 
k = 3.3, t = 0.01, and various |q|: (a) \a\ = 10 2 , (b) \a\ = 10 3 , 
(c) \a\ — 10 4 . The difference between exact and approximate 
results (both of them shown in the plot) is almost invisible. 
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FIG. 5. The quadrature uncertainty AX V versus |q| for 
k = 3.3, = 0, and various r: (a) r = 0.05, (b) r = 0.02, (c) 
t = 0.01, (d) r = 0.005. Exact results are plotted with solid 
lines and approximate ones with dashed lines. 



9 



